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ABSTRACT 



We present a method to calculate masses for components of both eclipsing and 
non-eclipsing binary systems as long as their apsidal motion rates are available. The 
method is based on the fact that the equation that gives the rate of apsidal motion 
is a supplementary equation that allows the computation of the masses of the com- 
ponents, if the radii and the internal structure constants of them can be obtained 
from theoretical models. For this reason the use of this equation makes the method 
presented here model dependent. 

We apply this method to calculate the mass of the components of the non-eclipsing 
massive binary system HD 93205 (03 V-l-08 V), which is suspected to be a very young 
system. To this end, we computed a grid of evolutionary models covering the mass 
range of interest, and taking the mass of the primary (Mi) as the only independent 
variable, we solve the equation of apsidal motion for Mi as a function of the age of 
the system. The mass of the primary wc find ranges from Mi = 60 ± 19 M© for ZAMS 
models, which sets an upper limit for Mi, down to Mi = 40 ± QM© for an age of 2 
Myr. Accordingly, the upper limit derived for the mass of the secondary (M2 = QMi) 
is M2 = 25 M0 is in very good agreement with the masses derived for other 08 V 
stars occurring in eclipsing binaries. 

Key words: binaries: eclipsing - stars: early type - stars: evolution - stars: funda- 
mental parameters - stars: individual (HD 93205) - stars: interiors 



1 INTRODUCTION 

■ The motion of the apse of a binary is mainly a direct con- 
sequence of the finite size of its components. If both stars 
were spherical objects and General Relativity corrections 
were negligible, they would move on a Keplerian, fixed or- 
bit. However, the presence of the companion object, and 
also its rotation, makes the structure of each star depart 
from spheres. In such a situation, there appears a finite 
quadrupolar (and higher) momentum to the gravitational 
field of each object that forces the orbit to modify the po- 
sition of the apse. This is an effect well known from long 
time ago (Cowling 1938; Sterne 1939). The rate of motion 
of the apse is dependent on the internal structure of each 
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component; thus, if we are able to determine the main char- 
acteristics of a binary, it provides an observational test of 
the theory of stellar structure and evolution (Schwarzschild 
1958; Kopal 1959). 

In spite of the age of the idea, apsidal motion of binary 
systems has been systematically studied only recently in a 
series of papers by Claret & Gimenez (e.g. Claret & Gimenez 
1993) and Claret (e.g. Claret 1995, 1997, 1998, 1999). Per- 
haps, one of the main reasons for such a situation is that 
the rate of motion of the apse is very dependent on the stel- 
lar structure. Thus, the apsidal motion test has been useful 
only recently because of the availability of accurate stellar 
evolutionary models such as those of Claret (1995). These 
authors have performed detailed stellar models computing 
the coefficients that determine the rate of motion of the apse 
and applied them to compare with observational data from 
eclipsing binaries. 

The detection of apsidal motion in non-eclipsing bi- 
nary systems is an elusive subject. It has to be determined 
through the time variation of the shape of the radial ve- 
locity curve due to change in the longitude of periastron. 
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Generally, the observed radial velocities in binary systems 
have large uncertainties that mask this effect in many cases. 
Moreover, the fact that times of light minima are usually 
determined with high precision for eclipsing binaries, acts 
as a selection effect in favour of detecting the motion of the 
line of apsides in such systems: a hundred cycles are usu- 
ally enough to notice the change in relative position of the 
secondary minimum with respect to the primary minimum. 
Observations over much longer periods of time are needed 
to find evidences of apsidal motion in systems where eclipses 
are not seen and the only observable effect is the change in 
shape of the radial velocity orbit. 

Empirical determinations of masses are scarce for early 
O-type stars, (e.g. Burkholder et al. 1997; Schonberner & 
Harmanec 1995). Few early O-type stars are known to be 
members of double-lined binaries, and from them, those 
showing eclipses or some kind of light variations that en- 
able the estimate of the orbital inclination (and then, abso- 
lute masses), are rare. The so-called "Mass Discrepancy", 
first described by Herrero et al. (1992) and recently re- 
viewed (Herrero, Puis, & Villamariz 2000), relates to the 
difference between the masses derived via numerical evo- 
lutionary models and those obtained from spectral analy- 
sis (plus model atmospheres) or binary star studies. This 
discrepancy, amounting to 50% in 1992, has been partially 
solved with the use of new evolutionary models that con- 
sider the effect of stellar rotation (Meynet & Maeder 2000) 
and new model atmospheres. But large differences between 
'predicted' and 'observed' masses are still present for the 
hottest and youngest (non evolved) stars. A recent study 
of the massive double-lined O-type binary system IfD 93205 
(Morrell et al. 2001) yields minimum masses of 31.5±1.1 Mq 
and 13. 3± 1.1 Mq for the 03 V primary and 08 V secondary 
components, respectively. This leads to a probable mass of 
~52-6OM0 for the 03 V star, if a mass value according to 
those derived for other 08V stars in eclipsing binaries is 
assumed for the secondary 08 V. This is much less than the 
80 to 100 Mq predicted from the position of this star on 
a theoretical HRD compared to stellar evolutionary tracks 
(see Fig. 7 in Morrell et al. 2001). 

Notably, HD 93205 is the first early O-type non- 
eclipsing ^ massive binary for which the rate of motion of the 
apse has been determined with some accuracy. Morrell et al. 
(2001) derived an apsidal motion period of 185 ± 16 years, 
that considering the orbital period of 6.0803 ± 0.0004 days, 
yields an apsidal motion rate rb — 0?0324 ± 0?0031 per or- 
bital cycle. This is very interesting because, from a mathe- 
matical point of view, the rate of motion of the apse provides 
another equation to be applied to the system apart from the 
standard ones. HD 93205 is an early-type massive short pe- 
riod binary in a highly eccentric orbit (e = 0.370 ± 0.005; 
Morrell et al. 2001), lying in the Carina Nebula, a galactic 
massive star forming region. Thus, we can assume its compo- 
nents are on, or very close to, the Zero Age Main-Sequence 
(ZAMS). Consequently, if we consider the evolutionary stage 
of HD 93205 as known, the equation of apsidal motion can 
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be written in a way that we can solve it for the mass of the 
primary star. The aim of the present paper is to detail this 
method and to apply it to HD 93205. 

The paper is organized as follows: in Section ^ we 
present the method we use to obtain the masses for compo- 
nents of non-eclipsing binary systems and we describe our 
evolutionary code and the calculations we have carried out. 
In Section we present a test of our method by applying 
it to some eclipsing binary systems. Section ^ is devoted to 
showing the results obtained for the massive binary system 
HD 93205 and finally, in Section |^ we give some concluding 
remarks about the implications of our results. 



2 COMPUTATIONAL DETAILS 

Here we present an original method, as far as the authors are 
aware, to calculate the masses of the components of binary 
systems provided the knowledge of the rate of motion of the 
apse and the evolutionary status of the stars. In addition, 
we describe the main characteristics of our evolutionary code 
and the calculation of the internal structure constants (ISC) 
on which the rate of the apsidal motion depends. 



2.1 Equations of apsidal motion and description 
of the method 

Sterne (1939) has shown that if the classical gravitational 
potential of each component of a binary system is expanded 
in a series of spherical harmonics, and terms up to the 
quadrupolar contribution are kept, then the rate of motion 
of the apse is given by |^ 



(1) 



where ri72 is the rate of the secular motion of the apse cal- 
culated considering only the quadrupolar contribution (the 
lowest order) of the gravitational potential; fc^j are the ISCs 
that depend on the int ern al mass distribution of the stars 
(see below. Subsection 2.2 for more details); i denotes the 
i-th multipolar momentum considered (i.e. i = 2 throughout 
this paper); whereas j denotes the component of the binary 
system. G is the gravitational constant, is the mean or- 
bital angular velocity, A is the semiaxis of the relative orbit, 
a\ and 02 are the mean radii of the stars, u)\ and UJ2 are 
their angular velocities of rotation. Mi and M2 are the stel- 
lar masses; and finally /2(e) and (72(e) are functions of the 
orbital eccentricity e given by 
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92 



(2) 
(3) 



In the following we shall reformulate Eq. (y) in order to 
write it to become an implicit equation for the mass of the 



^ It is assumed that rotation of both components is perpendicu- 
lar to the orbital plane. 
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primary star Mi. The semiaxis A is not directly known from 
observations, however the projected semiaxis D given by 



D — Asini 



(4) 



where sini is the sine of the inclination of the orbit, can be 
observationally assessed. Let us define the mass ratio Q and 
the angular velocities ratio qu: as 

In order to eliminate sini we can use the mass function, 
defined as 



/ 



Ml sin^ i 



Ml sin"' i 



(M1+M2)" {l + QY 



(6) 



which can be determined from observations (Batten 1973). 
In the same trend, projected tangential velocities vi — 
Vi sin i and V2 ~ V2 sin i are also observable quantities and 
their ratio q = V2/V1 can be used to eliminate q^j. 

One important point is that rotation modifies the inter- 
nal structure of the stars. In a recent paper. Claret (1999) 
has shown that, within the quasi-spherical approximation, 
rotation can be taken into account in the apsidal motion 
analysis simply by reducing the ISC ^2,1 by 



logfc2,i = log [k2,i]sph — 0.87Ai. 



(7) 



Here, [k2,i]sph denotes the ISC obtained from spherical mod- 
els and the parameter A; is defined by 



X^ = 



(8) 



where i denotes the component of the binary system and Vi, 
tti and Qi are, respectively, the tangential velocity, the radius 
and the surface gravity of the component. 

Up to this point, we have only considered the contribu- 
tions to the motion of the apse due to Newtonian gravity. 
However, it is well known that General Relativity predicts 
a secular motion of the apse which is independent of the 
classical contributions. The angular velocity of the apse due 
to General Relativistic effects ujcr is given by (Levi-Civita 
1937) 



^««= 6.36x10- ^^^ + ^^ 
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and incorporating both, rotation and relativistic effects we 
get after some algebraic manipulation 
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= 



where k2,i are the ISCs corrected by the effects of rotation. 
This is the fundamental equation for our purposes. 

As mentioned above, in Eq. ( |ll| ) some quantities are 
determined observationally (SI, ri72, e, D, Q, ui, /, g). On the 



other hand, k2,i and A;2,2 must be computed from evolution- 
ary models and, if we are dealing with non-eclipsing systems, 
as is the case for HD 93205, then the radii ai and 02 must be 
obtained from theoretical models as well. Now, if we assume 
that both components of the binary system have the same 
age, and we use the observational constraint Q = M2/M1, 
then A:2,i, fc2,2, Hi and 02 can be derived from evolutionary 
calculations as a function of Mi and the age of the system. 

The method presented here to calculate Mi is as fol- 
lows: we compute a grid of evolutionary models covering 
the range of masses of interest with a small mass step. Us- 
ing this grid, we construct isochrones starting at the ZAMS 
with a given time step and for each isochrone we seek the 
solution of Eq. (|ll|). In this procedure, the only independent 
quantity is Mi so when the solution of Eq. (hll) is found, the 
corresponding value of Mi is the mass of the primary star 
that corresponds to the age of the isochrone. 

Thus, for a given age of the system we have one solution 
for the mass of the primary star Mi, and using Q = M2/M1 
we can derive the mass of the secondary. However, the age of 
the system must be constrained by other means. In addition, 
the masses of the components that can be found with the 
present method are model dependent. Notice especially the 
sensitivity of the tidal and rotational terms in Eq. ( [ll| ) to the 
value of the stellar radius. It is clear that we need accurate 
stellar models in order to get a physically reliable value of 
Ml. 

Finally, it is worth mentioning that the solution of 
Eq. ( |ll| ) is subject to some constraints. Two of them were 
already mentioned: the age of both components in a binary 
system must be the same, and M2 = Q Mi. In addition, the 
value of the mass function imposes a minimum value for Mi 

[Ml]^^„ = fi^ + Qf■ (12) 

and also, the condition for the system to be detached 
ai+a2<D, (13) 
must be fulfilled. 



2.2 Evolutionary models and calculation of the 
ISCs 

As stated above, in order to solve Eq. (pT[), the ISCs and 
radii of both components must be obtained from evolution- 
ary calculations. This leads us to the necessity of having a 
set of evolutionary tracks of objects covering the range of 
masses expected for the components of the system. In the 
case of HD 93205, the primary 03 V star is candidate to be 
one of the most massive stars known. Thus, we have car- 
ried out calculations up to quite large stellar masses such as 
106 M0 because, as far as we are aware, there is no com- 
putation of the ISCs for such massive stars available in the 
literature. 

The calculations have been carried out with the stellar 
evolution code developed at La Plata Observatory. It is es- 
sentially the same code employed for studying white dwarf 
(see, e.g., Benvenuto & Althaus 1998) and intermediate mass 
stars (Brunini & Benvenuto 1997) and has been adapted for 
properly handling the case of massive stars. 

Let us briefiy describe the main ingredients of the code. 
The equation of state employed is that of OPAL (Rogers, 
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Swenson & Iglesias 1996). Radiative opacities are the lat- 
est version ol OPAL (Iglesias & Rogers 1996) while for low 
temperatures they are complemented with the Alexander & 
Ferguson (1994) molecular opacities. Conductive opacities 
and neutrino emission rates are the same as in Benvenuto 
& Althaus (1998). Nuclear reaction rates are taken from 
Caughlan & Fowler (1988) and weak electron screening is 
taken from Graboske et al. (1973). 

As we are dealing with massive stars, it is important to 
mention that we have accounted for the occurrence of over- 
shooting by employing the formalism described in Maeder & 
Meynet (1989). We have adopted the distance of overshoot- 
ing dov to be a fraction of the pressure scale height Hp at the 
canonical border of the convective zone: dov ~ oiovHp. Also, 
we allowed for mass loss following De Jager et al. (1986). 

Using the evolutionary code just described, we have cal- 
culated a set of evolutionary sequences covering the mass 
range from 4M0 to IO6M0 with a mass step of « 5%. We 
followed the evolution starting at the ZAMS till the deple- 
tion of hydrogen at the centre of the star. The initial helium 
content of our models is y = 0.275 and the adopted value 
for metallicity is Z = 0.02 while two values for overshooting, 
Uov = 0.25 and 0.40, were considered. 

After convergence of each model was reached, we solved 
the Clairaut-Radau differential equation (Sterne 1939) that 
accounts for the apsidal motion to the lowest (second) order 
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subject to the boundary condition 772 = at a = 0. In this 
expression a is the mean radius of a given equipotential, p(a) 
is the density at a and p(a) is the mean density interior to 
a, r]i is given by 

and the radius r of the distorted configuration and a are 
related by. 



(16) 



where Yi{a,9) describe the amplitude of the distortions. In 
order to integrate Eq. (^) we recall that p(a) and p{a) are 
provided by the structure of the evolving model. Close to 
the centre p{a), and consequently p(a) and 172(01), are ex- 
panded following an analogous treatment to that presented 
in Brooker & OUe (1955). The integration of Eq.(|l^) is 
started at the mesh point adjacent to the centre. Numer- 
ical integration is carried out with a standard Runge-Kutta 
routine (Press et al. 1986) up to the surface of the stellar 
model in order to get ?72 (oi). Then, the ISC k2,i is finally 
given by 



k2,i = 



3 - 772(a) 

4 + 2r;2(a) 



(17) 



3 A TEST OF THE METHOD EMPLOYING 
ECLIPSING BINARY SYSTEMS 

The method we are presenting here is, to our knowledge, 
original. In view of this fact, we have applied this method 



to some previously studied massive eclipsing binary systems 
with the aim of testing the method before applying it to HD 
93205. We have focused our attention on detached systems 
in which both components are massive stars in the main 
sequence (MS) and have a rather well measured apsidal mo- 
tion. We have finally selected the following systems: EM Car, 
QX Car, GL Car, Y Cyg and V478 Cyg. Observational pa- 
rameters for these systems are summarized in Table 0. In 
order to test the method, we compare the masses it yields 
with the observed ones, i.e. those obtained from the simul- 
taneous analysis of light and radial velocity curves of the 
systems. 

Components of a binary system must have the same age, 
so the first test we apply to the evolutionary code is that 
for each system considered there must be a single isochrone 
fitting the mass and radius of both components on the M — R 
plane. The results of our calculations for the choice aov ~ 
0.25 can be appreciated in Fig. |l| in which we show the 
mass and radius of the components of the selected binary 
systems. Note that for each system there is one isochrone 
that fits well both components, so the constraint that the 
ages of the components impose on evolutionary calculations 
is clearly satisfied by our models. In Fig. ^ we show the 
effective temperatures derived from our models for each star 
as a function of the observed effective temperature. Again, a 
good agreement between our evolutionary calculations and 
observations is found. As we have already stated, we have 
also considered a larger amount of overshooting, by fixing 
Uov to 0.40. We find no significant differences with the case 
of a smaller amount of overshooting, so we adopt the lower 
value (uov = 0.25) as the standard one in our calculations. 

Based on these preliminary results, we are confident 
that our models are appropriate for studying the evolution of 
stars in detached binary systems. Thus, we apply our mod- 
els to the study of apsidal motion through the calculation 
of the ISCs. Let us consider again Eq. (|l|). We can rewrite 
it in order to define a mean observational ISC k2.oba 



ZU2 -r , s 

— — K2,lC2,l + fc2,2C2,2 = K2,obs(C2,l + £2,2), 

and a mean theoretical ISC by 

Y _ fc2,lC2,l + fc2.2C2,2 

fe2,theo — \ . 

C2,l + C2,2 

where C2,i is given by 



C2,i 



15^^ /2(e) + -irn^92{e) 



M,G- 



(18) 
(19) 

(20) 



Even when the quotient VJ2 can be assessed from ob- 
servation, it is not possible to separate the contribution of 
each component to the rate of apsidal motion. Instead, we 
can determine k2,obs and it is this value that is currently used 
to contrast evolutionary models with observation. In Fig. ^ 
we show the theoretical values for k2 derived from our mod- 
els against the observed ones. We find that our models pre- 
dict mean ISCs that are in reasonable good agreement with 
the observed ones for the less concentrated models (those 
with a higher value of ^2). The most discrepant case we find 
is that of EM Car, which is the most evolved system con- 
sidered by us, as we find that its primary star has spent 
« 60% of its life on the MS. In this case, our models re- 
sult less concentrated than what we should expect from the 
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observed value k2,obs- However, as stated by Andersen & 
Clausen (1989), the apsidal motion rate for this system is 
based on observations covering only about 1 /6 of the apsidal 
motion period, thus the accuracy of apsidal motion parame- 
ters is still limited. In addition, information on its chemical 
composition is also missing. The theoretical values fc2,theo 
we find for the other systems are in good agreement with 
the observed values and, for QX Car and Y Cyg, also with 
those derived theoretically by Claret (1997). 

We present below some of the results of applying our 
method to the solution of Eq. (|ll]). First of all, let us em- 
phasize that for an assumed age of the system, the solution 
of Eq. ( [ll| ) is very well determined, i.e. only one solution is 
found as T2 is a very well-behaved, monotonously decreas- 
ing function of the independent quantity Mi. We illustrate 
this general behaviour with one example: in Fig. ^ we show 
T2 as a function of Mi for the case of V478 Cyg assuming an 
age of 6 Myr for both components. In view of these results, 
we find the method to be very reliable, from a mathematical 
point of view, in yielding a well determined value of M\. In 
Fig. P it is shown the mass Mi of the primary component of 
EM Car obtained as a function of the age of the system. The 
figure corresponds to the choice aov = 0.25 for overshooting. 
We find a good agreement between our theoretical predic- 
tion for Ml and its observed value for the whole range of ages 
considered, within a ilcr error. The upper limit for the age 
considered is given by the fact that, for larger ages. Mi falls 
below the minimum mass derived for this system from its 
mass function /. Results are very similar if an overshooting 
amount of aov = 0.40 is considered, the only main difference 
being that solution curve is slightly shifted to larger ages (« 
10%). Fig. ^ shows the results obtained for V478 Cyg again 
for aov ~ 0.25. For this binary system, an excellent agree- 
ment is achieved between our method and the observed mass. 
The same trend as before is found, i.e. the larger overshoot- 
ing shifts the solution toward ages about 10% larger. Note 
also that the larger the age considered, the smaller the mass 
of the primary (and also the mass of the secondary) that 
can account for the observed rate of apsidal motion. 



4 CALCULATION OF THE MASSES OF 

HD 93205 AND RELATED PARAMETERS 

From the results of the previous sections, we judge our 
method to be good enough to be employed in the mass es- 
timation of the components of HD 93205. As stated before 
in the Introduction, HD 93205 is a highly eccentric system, 
which strongly suggests that it must be very young. How- 
ever, the age estimates for such early O-type stars are very 
uncertain, either one tries to derive them considering the 
region in which the star is located, or comparing the star's 
position on the theoretical H-R diagram with isochrones 
calculated from evolutionary stellar models. HD 93205 be- 
longs to the open cluster Trumpler 16, the most massive 
stars of which have an age spread between 1 Myr and 2 Myr 
(DeGioia- Eastwood et al. 2001). Besides, there is evidence of 
ongoing star formation in the molecular cloud complex asso- 
ciated with the Carina Nebula (Megeath et al. 1996). Con- 
sequently, a lower limit to the age of the members of Tr 16 
cannot be established. On the other hand, de Koter, Heap, 



& Hubeny (1998) showed that if we increase by about 10% 
the effective temperature of 03-type stars, the age would 
decrease from 2 Myr to 1 Myr. Regarding the interpretation 
of theoretical isochrones for the most massive stars, these 
authors stated: "The derived Teff values are so similar be- 
cause the isochrone for ~ 2 Myr runs almost vertical and 
because the distance in temperature between the isochrones 
of 1 and 3 Myr is very small". 

Taking into account the problem in the age determi- 
nation described in the previous paragraph we choose to 
solve Eq. ( [ill ) ^'^^ ^ whole set of isochrones ranging from 
the ZAMS up to 2 Myr. We consider as zero age isochrone 
the one corresponding to the time when the stellar radius 
reaches its minimum value. In our models, this happens for 
ages of a few ten thousand years. 

In Fig.|^, we present the mass Mi of the primary compo- 
nent of HD 93205 as a function of the age of the system. Two 
curves are shown, each of them corresponding to a particu- 
lar choice of the overshooting parameter {uov = 0.25, 0.40). 
Let us emphasize that the amount of overshooting that ac- 
tually occurs is a rather uncertain quantity so we consider 
it as a free parameter and study its influence on the solu- 
tion of Eq. ([ll|). As can be seen from Fig. ^, for a given 
age. Ml is almost insensitive to our different choices of aov 
Both curves are almost overlapped over the whole range of 
ages, though differences tend to increase with age. This is 
not surprising because for a given mass the initial model (a 
ZAMS model) is the same in both cases so no initial discrep- 
ancy exists between them. As models evolve both sequences 
depart from one another and different internal mass concen- 
trations slowly arise. In view of this insensitivity, we shall 
concentrate ourselves on the case aov ~ 0.25 but there is no 
particular reason to prefer this value instead of the higher 
one. 

Let us consider again Fig. |^. The mass of the primary 
is a decreasing function of the age of the system. We find 
that its maximum value, corresponding to ZAMS models, 
is Ml = 60 ± 19M0. This is the upper limit for the mass 
All of the 03 V component of HD 93205. At increasing ages, 
it rapidly decreases and reaches 53 M0 at just 0.3 Myr and 
46.5M0 at 1 Myr approximately and finally Mi = 40±9 Mq 
at 2 Myr. Within observable quantities, the main source of 
uncertainty in determining Mi is the apsidal motion rate 
(known up to a 9% accuracy) and to a smaller extent the 
projected semi-axis and the projected rotational velocities, 
so better determinations of these quantities (especially the 
apsidal motion rate) are needed in order to decrease the error 
in the determination of Mi. We recall here that the apsidal 
motion rate is a critical parameter because the necessity of 
very long time baseline (decades) of high-quality observa- 
tions. Having the mass Mi determined, it is straightforward 
to calculate the mass M2 of the secondary if we recall (Table 
^ see Morrell et al. 2001 for further details) that the mass 
ratio Q = M2/M1 for HD 93205 is 0.423 ± 0.009. We find 
that M2 ranges from 25.3 ± 8 Mq at the ZAMS down to 
17 ± 4M© if a rather large value of 2 Myr is adopted for the 
age of the system. These mass values are in good agreement 
with those expected for an 08 V star like this one (consider, 
particularly, the well known short-period eclipsing binary 
EM Car, whose primary component is an 08 V and its mass 
is 22.89 ± O.32M0, Andersen & Clausen 1989). Once Mi is 
determined it is easy to obtain the inclination i of the orbit 
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Table 1. Astrophysical parameters for selected test systems 
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from Eq. In Fig. |^ it is shown the resulting inclination 
from the set of calculations corresponding to ao„ = 0.25. It 
can be seen that the inclination of the system increases as 
age does. This is a direct consequence of the behaviour of 
Ml (which decreases as age increases) but it is worth noting 
that within the whole range of ages considered the resulting 
inclination does not allow eclipses to occur. Indeed, if we 
assume that HD 93205 is not older than 2 Myr we find that 
54° < i < 68° , in coincidence with Antokhina et al. (2000) 
who found a most probably value of i = 60° . 

However, a problem arises when we try to compare the 
luminosity derived from the corresponding models to the 
observed value for HD 93205. Let us explain this with an ex- 
ample: if we consider the 60 M© model, it predicts, for zero 
age, a radius Ri = 10.7 Rq, and a logTcff — 4.68, result- 
ing in a luminosity, logL = 5.72 Lq. This corresponds to a 
bolometric magnitude, Mboi = —9.55, which is almost one 
magnitude fainter than Mboi = —10.41 derived by Morrell 
et al. (2001) from the visual magnitude of the 03 V compo- 
nent of HD 93205, the distance modulus of 12.55 obtained 
by Massey & Johnson (1993) for Tr 16, and the bolometric 
correction (BC) for an 03 V star taken from the calibration 
by Vacca, Garmany & ShuU (1996). This large disagreement 
between the expected and observed bolometric magnitudes, 
points to a large error in some (or any) of the involved as- 
sumptions. If the distance modulus is right, then we can 
suspect that the BC must be wrong by about one magni- 
tude. On the other hand, the distance modulus of the Carina 
Nebula is still a matter of discussion. Distance modulus of 
the order of that derived by Massey & Johnson (1993) arise 
from the consideration of color-magnitude diagrams for the 
stellar component of the clusters. Some other independent 
determinations, like the recently obtained by Davidson et 
al. (2001) from kinematic study of the Homunculus nebula 



Table 2. Observed parameters for HD 93205 



ai sinj [km] 
a2 sini [km] 
Ki [km s-i] 
K2 [km s-i] 
P [days] 
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Ml sin^ i [Mq' 
M2 sin^ « [Mq' 
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m [° days-l] 
Vi sini [km s~ 
V2 sini [km s~ 



(1.03 ±0.02) 10'' 
(2.44 ±0.02) 10'' 
132.6 ±2.0 
313.6 ± 1.8 
6.0803 ± 0.0004 
0.370 ± 0.005 
31.5 ± 1.1 
13.3 ± 1.1 
0.423 ± 0.009 
0.00533 ± 0.00051 
135 
145: 



Morrell et al. (2001) 



Howarth et al. (1997) 



surrounding Eta Car, give distance modulus as low as 11.76, 
which would significantly decrease the referred discrepancy. 
But, if we suppose this last distance modulus to be correct, 
then all of the stars in Tr 16 will have Mv about 0.8 mag- 
nitudes fainter than the values accepted to date. Here we 
arrive at a point whose importance is obvious for many as- 
trophysical issues, and deserves to be carefully studied. The 
referred discrepancies might also arise in a combination of 
different sources of error (BCs, distances, and adopted abso- 
lute magnitude scale for ZAMS stars) . A detailed discussion 
of these issues will be presented in a forthcoming paper. 

Finally, let us comment briefly that tidal contribution to 
the apse motion of HD 93205 is the most important, rang- 
ing from about 60% at the ZAMS to 70% at 2 Myr. The 
rotational contribution ranges from 30% to 20% and the 
relativistic one is almost constant and approximately 10% 
of the apsidal motion rate is due to this effect. In this sense, 
HD 93205 could be classifled as a relativistic binary system 
(Claret 1997). 
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5 CONCLUSIONS 

We present a method to calculate masses for components of 
non-eclipsing binary systems if their apsidal motion rate is 
provided. The method consists in solving Eq. ( pj| ) if the ra- 
dius and the internal structure constant of each component 
can be obtained from a grid of stellar evolution calculations. 
In order to test this method, we have selected some eclipsing 
binary systems and have derived the masses of their compo- 
nents. A very good agreement was achieved between masses 
obtained with our method and those derived from the anal- 
ysis of their radial velocity and light curves. 

The main goal of this article, besides presenting the 
method, is to calculate the masses of the components of 
HD 93205. This is an 03V -I- 08V system. Its 03V com- 
ponent has the earliest known spectral type of a normal 
star found in a double-lined close binary system, thus poten- 
tially being a very massive star. Although HD 93205 is not 
an eclipsing binary, Morrell et al. (2001) have measured its 
apsidal motion rate and found it to be ri? = 0?0324±0?0031 
per orbital cycle so we have been able to apply the method 
presented here to this system. The resulting mass of the pri- 
mary star (Ml) is obtained as a function of the assumed 
age of the system. HD 93205 is a highly eccentric system 
(e = 0.370 ± 0.005) which suggests a very low age. How- 
ever, we do not adopt a particular value for the age as its 
determination is quite uncertain, and prefer to consider a 
range of ages starting at the ZAMS. We find that for zero 
age models the resulting mass is Mi = 60 ± 19 Mq and 
that it monotonously decreases as age is increased (Fig. |^, 
reaching Mi = 40 ± 9 M© at 2 Myr. Now, if we take into 
account the mass ratio Q = 0.423 for HD 93025, the mass 
of the secondary lies in the range M2 = 25.3 — 17 M© for 
this range of ages. It is worth mentioning again that these 
M2 values are in good agreement with the masses derived 
for other 08 V stars in eclipsing binaries such as the well 
studied system EM Car (Andersen & Clausen 1989). The 
mass value derived for Mi is also in the range (52 - 60 M0) 
obtained from the observed Q assuming a "normal" mass 
for the 08V secondary component (i.e. 22 - 25 M0). In 
addition, we have estimated the inclination of the system 
through Eq. ^ and the results obtained (Fig. ^) are con- 
sistent with the non-eclipsing condition of HD 93205. 

Our results corresponding to zero age give an upper 
limit to the mass of the 03V component of HD 93205, a 
result that places a strong constraint to the masses of the- 
oretical stellar models for the most massive stars. Also, the 
luminosity derived from the stellar models for the 03V 
component rises a problem when compared with the ob- 
served value, being the theoretical Mboi almost one magni- 
tude fainter than the value derived from the observations. 
This discrepancy raises the need of reviewing both the dis- 
tance and BC scales for the earliest type ZAMS stars, a 
subject that will be addressed in the near future. 
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Figure 1. The radius vs. mass relationstiip for tiie components of binary systems EM Car, QX Car, Gl Car, Y Cyg and V478 Cyg 
together with their corresponding error bars. Solid lines represent our theoretical isochrones for 7.8, 6.3, 4.3 and 1.8 Myr (from right to 
left). For each system considered there is one isochronc that fits both components. 
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Log 



Figure 3. Comparison between theoretical and observed k2 for the components of the same binary systems included in Fig. |j 1 o" error 
bars for k2^oba are also shown. Notice that theoretical models are slightly less concentrated than indicated by observations. 
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Figure 4. T2 as a function of Mi for tiie case of V478 Cyg, assuming an age of 6.0 Myr for the system and setting aov = 0.25. 
that the solution value M\ for which T2 = is very well determined, since T2 is a monotonously decreasing function of M\. 
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Figure 5. The mass of the primary star of EM Car system assuming oiov = 0.25 as a function of its age. Horizontal solid line 
corresponds to the preferred value deduced from radial velocity and light curves. Short dashed horizontal lines represent the uncertainty 
in this observational value. The other solid and short dashed lines represent the preferred value and its Icr uncertainty respectively, 
deduced from the apsidal motion rate. 
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Figure 6. Same as Fig. |^ but for V478 Cyg 
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Figure 7. The mass of the components of HD 93205 deduced from its apsidal motion rate as a function of its assumed age. In each 
figure, sohd lower (upper) line corresponds to the preferred value assuming Oov = 0.25 (ckou = 0.4). Short dashed lines represent its Icr 
uncertainty. For more details, see text. 
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Figure 8. Inclination of HD 93205 as a function of its age. For more details, see text. 



© 2001 RAS, MNRAS 000, 0-| 



